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A quantum mechanical bound for CHSH-type 
Bell inequalities 
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Abstract Many typical Bell experiments can be described as follows. A source 
repeatedly distributes particles among two spacelike separated observers. Each of 
them makes a measurement, using an observable randomly chosen out of several 
possible ones, leading to one of two possible outcomes. After collecting a sufficient 
amount of data one calculates the value of a so-called Bell expression. An important 
question in this context is whether the result is compatible with bounds based on 
the assumptions of locality, realism and freedom of choice. Here we are interested 
in bounds on the obtained value derived from quantum theory, so-called Tsirelson 
bounds. We describe a simple Tsirelson bound, which is based on a singular value 
decomposition. This mathematical result leads to some physical insights. In par¬ 
ticular the optimal observables can be obtained. Furthermore statements about the 
dimension of the underlying Hilbert space are possible. Finally, Bell inequalities 
can be modified to match rotated measurement settings, e.g. if the two parties do 
not share a common reference frame. 


1 Introduction 

Since the advent of quantum theory physicists have been struggling for a deeper un¬ 
derstanding of its concepts and implications. One approach to this end is to carve out 
the differences between quantum theory and “classical” theories, i.e. to explicitly 
point to the conflicts between quantum theory and popular preconceptions, which 
evolved in each individual and the scientific community from decoherent macro¬ 
scopic experiences. Plain formulations of such discrepancies and convincing exper¬ 
imental demonstrations are crucial to internalizing quantum theory and replacing 
existing misconceptions. For this reason the double-slit-experiments (and similar 
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Fig. 1 Two parties, Alice (A) and Bob (B), perform a Bell experiment. Both of them receive parts 
of a quantum system from the source (S). They randomly choose a measurement setting, denoted 
by x= 1,2, andy = 1,2, and write down their outcomes a = — 1 or 1 and b = — 1 or 1, 

respectively. The experiment is repeated until the accumulated data is analyzed according to the 
text. Angles of 45 degrees in the space-time-diagram correspond to the speed of light. The future 
light cones of A, B and S show, that the setting choice and outcome of one party cannot influence 
the other and that A and B also cannot influence any event inside the source. 


experiments with optical gratings) CD 12 0110, which expose the role of state su¬ 
perpositions in quantum theory, are so very fascinating and famous. Other examples 
of “eye-openers” are demonstrations of tunneling (6] [7] PP- 33-12], the quantum 
Zeno effect HI and variations of the Elitzur-Vaidman-scheme HU HMD, to pick 
just a few. 

Bell experiments fT2lll3 , 14 f5|- which show entanglement in a particularly strik¬ 
ing way, belong to this list. Informally, entanglement is the fact that in quantum 
theory the state of a compound system (e.g. two particles) is not only a collection of 
the states of the subsystems. This fact can lead to strong correlations between mea¬ 
surements on different subsystems. Before going into more detail here, we would 
like to note that the described differences between the relatively new quantum theory 
and our old preconceptions are obvious starting points when to look for innovative 
technologies which were even unthinkable before. This is in fact a huge motivation 
for the field of quantum information, where Bell experiments play a central role. 


1.1 Bell experiments bring three fundamental common sense 
assumptions to a test 

The idea of Bell was to show that some common sense assumptions lead to predic¬ 
tions of experimental data which contradict the predictions of quantum theory. In the 
following we employ a black box approach to emphasize that this idea is completely 
independent of the physical realization of an experiment. For example the measure¬ 
ment apparatuses get some input (an integer number which will in the following be 
called “setting”) and produce some output (the “measurement outcomes”). We refer 
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readers preferring a more concrete notion to Section |~L2| where physical implemen¬ 
tations and concrete measurements are outlined. 

In the present paper we consider the following (typical) Bell experiment, see also 
Figure [T| There are three experimental sites, two of which we call the parties Al¬ 
ice (A) and Bob (B), and the third being a preparation site which we call source 
(S). Alice and Bob have a spatial separation large enough such that no signal can 
travel from one party to the other at the speed of light during the execution of our 
experiment. The source is separated such that no signal can travel from A or B to it 
at the speed of light before it finishes the state production. The importance of such 
separations will become clear later. 

The source produces a quantum system, and sends one part to Alice and one to Bob. 
We will exemplify this in Section [L2| A and B are in possession of measurement ap¬ 
paratuses with a predefined set of different settings. In each run they choose the set¬ 
ting randomly, e.g. they turn a knob located at the outside of the apparatus, measure 
the system received from the source and list the setting and outcome. In the present 
paper the measurements are two-valued and the outcomes are denoted by — 1 and 
+ 1. Let Mi and M 2 be the number of different measurement settings at site A and B, 
respectively. We label them by x = 1,2,3,..., Mi for Alice and y = 1,2,3,..., Mo for 
Bob. This preparation and measurement procedure of a quantum system is repeated 
until the amount of data suffices to estimate the expectation value of the measured 
observables, up to the statistical accuracy one aims at. The expectation value of an 
observable is the average of all possible outcomes, here ±1, weighted with the cor¬ 
responding probability to get this outcome. 

Let us sketch the preconceptions that are jointly in conflict with the quantum theo¬ 
retical predictions for Bell tests. These are mainly three concepts: Locality, realism 
and freedom of choice. This forces us to question at least one of these ideas, be¬ 
cause any interpretation of quantum theory, as well as any “postquantum” theory, 
cannot obey all of them. We invite the reader to pick one to abandon while reading 
the following descriptions. Do not be confused by our comparison with the textbook 
formalism of quantum theory: so far you are free to choose any of them. 

Locality is the assumption, that effects only have nearby direct causes, or the other 
way around: any action can only affect directly nearby objects. If some action here 
has an impact there, then something traveled from here to there. And, according to 
special relativity, the speed of this signal is at most the speed of light. In our setup, 
this means that whatever Alice does cannot have any observable effect at Bob’s site. 
In particular, the measurement outcome at one side cannot depend on the choice of 
measurement setting at the other site. While the formalism of quantum theory has 
some “nonlocal features”, e.g. a global state, it is strictly local in the above sense, 
because any local quantum operation on one subsystem does not change expectation 
values of local observables for a different subsystem. 

Realism is the concept of an objective world that exists independently of subjects 
(“observers”). A stronger form of realism is the “value-definiteness” assumption 
meaning that the properties of objects always have definite values, also if they are 
not measured or even unaccessible for any observer. It seems to be against common 
sense to assume that objects cease to have definite properties if we do not measure 
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them any longer. In particular the natural sciences were founded on the assumption, 
that nature and its properties exist independently of the scientist. In our setup realism 
implies, that the measurement outcomes of unperformed measurements (in uncho¬ 
sen settings) have some value. We do not know them, but we can safely assume that 
they exist, give them a name and use them as variables. If possible outcomes are 
— 1 and +1, for example, we might use that the outcome squared is 1 in any of our 
calculations. In general the (usual) formalism of quantum theory does not contain 
definite values for measurement outcomes independent of a measurement. 

Freedom of choice, which is also sometimes called the free will assumption, means 
it is possible to freely choose what experiment to perform and how. Because this 
idea is elusive, we are content with a decision that is statistically independent of any 
quantity which is subject of our experiment. The idea of fate seems to be tempting 
to many people. However, dropping freedom of choice makes science useless. Just 
imagine you “want” to investigate the question whether a bag contains black balls 
but your fate is to pick only white balls (and put them back afterwards), even though 
there are many black balls inside. In our setup, freedom of choice implies, that A’s 
and B’s choice of measurement setting does not depend on the other’s choice or the 
outcomes. In quantum theory, there is freedom of choice in the sense that random 
measurement outcomes of some other process can be used to make decisions. 

If you decided that you preferably take leave of locality you are in good company. 
Many scientists conclude from Bell’s theorem, that the locality assumption is not 
sustainable. This is particularly interesting when you consider the above compari¬ 
son with the standard textbook formalism of quantum theory, which is apparently 
not realistic but local in the described sense. The fact that in this context many scien¬ 
tists speak about “quantum nonlocality” thus leads to controversy Ifl6l . We therefore 
want to stress again, that the experimental contradiction only tells us that at least one 
of all the assumptions that lead to the predictions needs to be wrong. We cannot de¬ 
cide which assumption is wrong from Bell’s theorem alone. 

We now focus on a tool to show the contradiction in the described experiment be¬ 
tween the above assumptions and quantum theory, the so called Bell inequalities. 
These are inequalities of measurable quantities which are (mainly) derived from lo¬ 
cality, realism and freedom of choice and therefore hold for all theories which obey 
these principles, while they are violated by the predictions of quantum theory. We 
consider a special kind of Bell inequalities which are linear combinations of joint 
expectation values of Alice’s and Bob’s observables. The joint expectation value 
of the two observables of Alice and Bob is the expectation value of the product of 
the measurement outcomes, which again takes values ±1. It depends on the setting 
choice x at Alice’s site and y at Bob’s site and we denote it by E(x,y). If we denote 
the (real) coefficient in front of the expectation value E{x,y ) as g x>y , then we can 
write such Bell inequalities as 


M\ M 2 

E Y,8x, y E(x,y)<B g , 


x=ly=l 


(1) 
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where the bound B g depends only on the coefficients g XJ . These coefficients form a 
matrix g which has dimension M\ x M?. Any real matrix g defines a Bell inequality 
via Eq. (TTl). The may be most famous example is the Clauser-Horne-Shimony-Holt 
(CHSH) Q7) inequality, which reads 

£(1,1)+£(1,2)+£(2,1)—£(2,2) <2. (2) 

Here the corresponding matrix g is 

Due to its prominence we call the class of Bell inequalities in the form of Eq. G 
CHSH-type Bell inequalities. For completeness we sketch the derivation of B g . It 
turns out that it suffices to consider deterministic outcomes only, as a probabilistic 
theory, where the outcomes follow some probability distribution, cannot achieve a 
higher value in Eq. (JTJ: it can be described as a mixture of deterministic theories and 
the value of Eq. G is the sum of the values for the deterministic theories weighted 
with the corresponding probability in the mixture. For deterministic theories the 
expectation value is merely the product of the two (possibly unmeasured) outcomes 
a of Alice and b of Bob, which we are allowed to use when assuming realism. Due 
to locality a only depends on the setting x of Alice, which has no further dependence 
due to freedom of choice. Analogously b depends only on the setting y of Bob, which 
in turn has no further dependence. Thus the expectation value is 

E(x,y)=a(x)b(y). (4) 

Now we can calculate B g by maximizing Eq. G over all possible assignments of 
— 1 and +1 values to a(x') and b(y). In Eq. (|2j» the maximal value is B g = 2, which 
is achieved for a(l) = a{ 2) = b( 1) = b( 2) = 1, for example. Note that the sign of 
£(2,2) cannot be changed independently of the other three terms, because £(1,2) 
and £(2,1) contain b( 2) and a( 2), respectively. 

We point out that any function that maps the probabilities of different measurement 
outcomes to a real number may be used to derive Bell inequalities, and different 
types of Bell inequalities can be found in the literature (e.g. my However, here 
we focus on Bell inequalities of the form of Eq. G- 


1.2 The CHSH inequality can be violated in experiments with 
entangled photons 

We recapitulate some basics of quantum (information) theory. Analogously to a 
classical bit the quantum bit, or qubit, can be in two states 0 and 1, but additionally 
in every possible superposition of them. Mathematically this state is a unit vector in 
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the two-dimensional Hilbert space (a vector space with a scalar product) C 2 spanned 
by the basis vectors 

0:= (J) and 1 := (5) 

An example of a superposition of these basis states is yr = (0 +1). Any observ¬ 

able on a qubit with outcomes +1 and — 1 can be written as 



Ox <7y Or 


where the vector a = (a x ,a y ,a z ) T (here 1 denotes transposition) defines the mea¬ 
surement direction and the matrices cr T , oy and a z are called Pauli matrices. The ex¬ 
pectation value of this observable given any state y/ can be calculated as E = y/^Ay/ 
(here ' denotes the complex conjugated transpose), which is between —1 and +1. 
Any quantum mechanical system with (at least) two degrees of freedom can be used 
as a qubit. In the present context the spin of a spin-^-particle, two energy levels of 
an atom and the polarization of a photon are important examples of qubits. The spin 
measurement can be performed using a Stern-Gerlach-Apparatus fl9l . the energy 
level of an atom may be measured using resonant laser light, or the polarization of 
a photon can be measured using polarization filters or polarizing beam splitters. 

The Hilbert space of two qubits is constructed using the tensor product, i.e. C 2 0 
(D 2 = C 4 . The tensor product of two matrices (of which vectors are a special 
case) is formed by multiplying each component of the first matrix with the com¬ 
plete second matrix, such that a bigger matrix arises. The state of the compos¬ 
ite system of two qubits in states (p A = and (j) B = then reads 

<j) AB = (j) A ®(j) B = {(l)f(j) B 7 (l)f(j) B : (l) A (j) B ,(l) A (j) B ) T . The states of such composite sys¬ 
tems might be superposed, which leads to the notion of entanglement. 

Out of several physical implementations of the CHSH experiment we sketch the 
ones with polarization entangled photons (see 1201 1. We identify 0 with the horizon¬ 
tal and 1 with the vertical polarization of a photon. Nonlinear processes in special 
optical elements can be used to create two photons in the state 


<j>+ = -^(l,0,0,l) r , (7) 

i.e. an equal superposition of two horizontally polarized photons and two vertically 
polarized photons. The measurements of Alice and Bob in setting 1 and 2 are 

A\ = cos(2 x 22.5°)oy-f-sin(2 x 22.5 °)ct z , (8) 

A 2 = cos(—2 x 22.5°)oy-t-sin(—2 x 22.5°)a ; , (9) 

B\ = cos(2 x 0 °)cj v + sin(2 x 0°)oy (10) 

and Bi = cos(2 x 45°)oy-t-sin(2 x 45°)cr : , (11) 
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respectively. Here the angles are the angles of the polarizer and the factor 2 is due to 
the fact that in contrast to the Stern-Gerlach-Apparatus a rotation of the polarizer of 
180° corresponds to the same measurement again. One can now calculate the value 
of Eq. [2] 

£(1,1)+£(1,2)+£(2,1)-£(2,2) = (j>^ + (Ai<S>Bi)(j> + + 0+(A[ ®£ 2 )0+ 

+0 1 ((A2®£i)0 + -0+(A2(8>£ 2 )0+ 

_L _L 

V2 + V2 + V2 V V2 ) 

= 2V2. (12) 

The value 2\/2 « 2.82 is larger than 2 and therefore the CHSH inequality is violated. 
One can ask whether it is possible to achieve an even higher value, e.g. when using 
higher-dimensional systems than qubits, because at the first glance a value of up to 
four seems to be possible. This question is addressed in the following sections (the 
answer, which is negative, is given in Section [2+|. 


1.3 The quantum analog to classical bounds on Bell Inequalities 
are Tsirelson Bounds 

Analogously to the “classical” bound one can ask for bounds on the maximal 
value of a Bell inequality obtainable within quantum theory, so-called Tsirelson 
bounds ED, and the observables that should be measured to achieve this value. In 
other words: which observables are best suited to show the contradiction between 
quantum theory and the conjunction of the three discussed common sense assump¬ 
tions. This question, which is also of some importance for applications of Bell in¬ 
equalities, is the main subject of the present essay. 

The scientific literature contains several approaches to derive Tsirelson bounds, 
some of which we want to mention. The problem of finding the Tsirelson bound 
of Eq. 0 can be formulated as a semidefinite program. Semidefinite programming 
is a method to obtain the global optimum of functions, under the restriction that the 
variable is a positive semidefinite matrix (i.e. it has no negative eigenvalues). This 
implies that well developed (mostly numerical) methods can be applied ll22l [23l. 
The interested reader can find a Matlab code snippet to play around with in the ap¬ 
pendix [4] Furthermore there has been some effort to derive Tsirelson bounds from 
first principles, amongst them the non-signalling principle ll24ll . information causal¬ 
ity l25l and the exclusivity principle j26l , 

The non-signalling principle is satisfied by all theories, that do not allow for faster- 
than-light communication. Information causality is a generalization of the non¬ 
signalling principle, in which the amount of information one party can gain about 
data of another is restricted by the amount of (classical) communication between 
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them. The exclusivity principle states, that the probability to see one event out of a 
set of pairwise exclusive events cannot be larger than one. 


2 The singular value bound 

Here we will discuss a simple mathematical bound for the maximal quantum value 
of a CHSH-type Bell inequality defined via a matrix g , which we derived in ||27|| . 
While it is not as widely applicable as the semidefinite programming approach, it is 
an analytical expression which is easy to calculate and it already enables valuable 
insights. For “simple” Bell inequalities, like the CHSH inequality given above, it is 
sufficient to use the method of this paper. 

We will make use of singular value decompositions of real matrices, a standard tool 
of linear algebra, which we now shortly recapitulate. 


2.1 Any matrix can be written in a singular value decomposition 

A singular value decomposition is very similar to an eigenvalue decomposition, in 
fact the two concepts are strongly related. Any real matrix g of dimension M\ x M 2 
can be written as the product of three matrices V, S, W T , i.e. 

g = VSW T , (13) 

where these three matrices have special properties. The matrix V is orthogonal, i.e. 
its columns, which are called left singular vectors, are orthonormal. It has a dimen¬ 
sion of M\ x Mi. The matrix S is a diagonal matrix of dimension M 1 x M 2 , which 
is not necessarily a square matrix. Its diagonal entries are positive and have non¬ 
increasing order (from upper left to lower right). They are called singular values of 
g. The matrix W is again orthogonal. It has dimension Mi x M 2 and its columns are 
called right singular vectors. 

The largest singular value can appear several times on the diagonal of S. We call the 
number of appearances the degeneracy d of the maximal singular value. Due to the 
ordering of 5, these are the first d diagonal elements of S. Here we note the concept 
of a truncated singular value decomposition: instead of using the full decomposition 
one can approximate g by using only parts of the matrices corresponding to, e.g., 
the first d singular values (i.e. only the maximal ones). These are the first d left and 
right singular vectors, and the first part of S , which is just a d x d identity matrix 
multiplied by the largest singular value. Since these matrices play an important role 
in the following analysis we will give them special names: V^, and W- d> . All 
these matrices are depicted in Figure [2] 

The matrix g maps a vector v to a vector gv which, in general, has a differ¬ 
ent length than v. Here the length is measured by the (usual) Euclidean norm 
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Fig. 2 The matrices involved in the singular value decomposition of a general real M\ x Mi matrix 
g: V and W are orthogonal matrices, S is diagonal. V and W contain the left and right singular 
vectors, respectively, as columns, and S contains the singular values on its diagonal. The shaded 
parts belong to a truncated singular value decomposition of g. We denote the parts corresponding 
to the maximal singular value as S 1 -' 1 ' 1 and 


| |v| I 2 = y v i + v 2 + ••• + v m 2 - The possible stretching factor for all vectors 

v is a property of the matrix: its matrix norm induced by the Euclidean norm. The 
value of this matrix norm coincides with the maximal singular value .S’j 1 . We can 
therefore express the maximal singular value using 

5,1 = ma « 7ETT =: Hsll 2 ' (14) 

VG1R M 2 ||V||2 

The notation | |g| I 2 for the maximal singular value of g is more convenient than Si 1 , 
as it contains the matrix as an argument. 


2.2 The singular value bound is a simple Tsirelson bound 

It turns out that the matrix norm of g , i.e. its maximal singular value, leads to an up¬ 
per bound on the quantum value for a Bell inequality, defined by g via Eq. Q. This is 
the central insight of this essay. It is remarkable that a mathematical property, solely 
due to the rules of linear algebra, leads to a bound for a physical theory, here the 
theory of quantum mechanics. With the definition of the matrix norm given above, 
we can now write this singular value bound of g, a simple Tsirelson bound E71 , It 
reads 

M\ M 2 

Tj TjSxi^ 2 E{xi,x 2 ) < VMMM 2 , (15) 

x— 1 y— 1 

where E now denotes the expectation value of a quantum measurement in setting 
*1 and X 2 - Eq. is the central formula of this essay. Note that this bound is not 
always tight, i.e. there exist examples where the right hand side cannot be reached 
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within quantum mechanics. However for many examples it is tight. The proof of 
this bound is sketched in Appendix [3] 

We now calculate this bound for the CHSH inequality given in Eq. |[2}. We see, that 
here the matrix of coefficients is 


8 = 


1 

1 



i 

y/2 



V2 0 \ (l 0\ 

0 V 01 /' 

r w t 


(16) 


It is easy to check that the given decomposition of g is a singular value decomposi¬ 
tion, i.e. V, S and W have the properties described above. From this we read, that the 
maximal singular value of g is ||g ||2 = \/2. Then Eq. (151 tells us, that the maximal 
value of the CHSH inequality (Eq. Q) within quantum theory is not larger than 
2y/2, a value which can also be achieved when using appropriate measurements and 
states (see Section |~i~2]) 


2.3 Tightness of the bound can be checked efficiently 

We already mentioned that the inequality is not always tight, i.e. sometimes it is 
not possible to find observables and a quantum state such that there is equality. From 
the derivation of Eq. m sketched in Appendix [3] one understands, why this is the 
case. The value y/MiAfWsWi is achieved if and only if there exists a right singular 
vector v to the maximal singular value and and a corresponding left singular vector 
w which fulfill further normalization constraints. 

It is common to denote the element in the /-th row and /'-th column of a matrix A 
as A ij. We will extend this notation to denote the whole /-th row by A,-* and the 
whole j -th column by A*;, i.e. the * stands for “all”. For example, the /-th M\ +M 2 
dimensional canonical basis vector, with a one at position / and 0 everywhere else, 
can then be written as 

With this notation at hand we write down the normalization constraint from above 
as the system of equations 


T T/(^) 


and 


a‘V , 


Mi T (d) 
TT a Wy* 
M 1 7 


= 1 for x = 1 , 2 ,..., Mi 


= 1 fory = 1,2,...,M 2 , 


(17) 

(18) 


where the d x d! matrix a is the unknown. The bound in Eq. (15 1 is tight if and 
only if such matrix a solving this system of equations can be found. Here d is the 
degeneracy of the maximal singular value of g and d', the dimension of the vectors 


T ( d) t id) 

\ x = a Vy* and w v = a Vf , is a natural number. The steps leading to Eqs. (17 1 

and (18 1 can be found in the supplemental material of ll27l . Because Eqs. (| 1 1\ and 


















A quantum mechanical bound for CHSH-type Bell inequalities 


11 


( [T8| are quadratic in a it may not be obvious how to solve it. In ll27l we described an 
algorithm to solve the above system of equations in polynomial time with respect to 
the size of g. The interested reader may also find a Matlab snippet in the Appendix[4] 
Often the solution a is obvious, e.g. when it is proportional to the identity matrix. 


2.4 Optimal measurements are obtained from the SVD 


From the previous considerations we understand that the existence of the unit vec- 

T (d) t (d) 

tors v A - = a V„ and w y = a Vy* , i.e. the existence of the matrix a that allows 
this normalization, is crucial to the satisfiability of the singular value bound. Fur¬ 
thermore they have a physical meaning, because they are related to the observables 
in the following way. 

Let us again consider the example of Eq. (|2j, with the singular value decomposition 


g = 


1 1 
1 -1 



(19) 


which we repeat from Eq. ( p~6| ). The multiplicity d = 2 of the maximal singular value 
\fl equals the number of measurement settings M\ and M 2 , so each of the rows vjf ! 
and wjf are already normalized due to orthogonality of V and W. Therefore we 


can choose a = to solve Eqs. (17 1 and (181. We then have vi = (1, \) T /\/2, 
\2 = (1, — \) T /\/2, Wi = (l,0) r and W 2 = (0, lj' . We are looking for a state and ob¬ 
servables such that E(x,y) = \ x ■ w v , which is always possible to find (see Tsirelson’s 
theorem, Appendix[3]). 

Consider for example two spin-j particles in the state (f> + = ^(1,0,0, l) r . Alice 
and Bob can measure their particles’ spin with Stern-Gerlach apparatuses along any 
orientation in the x-z-plane. The observable of Alice corresponding to a measure¬ 
ment along the direction (a x ,a z ) T is 


A — a x 




( 20 ) 


where the matrices are two of the so-called Pauli matrices. Bob’s measurement reads 
analogously. The reader can easily verify that the expectation value of the joint 
observable A (g> B is given by 

<t>\(A®B)<i> + = a-b. (21) 

Therefore optimal measurement directions leading to equality in Ineq. ( fj~5] l are given 
by \ x and w v . For this reason we will call \ x and w v the measurement directions, 
even though they can have a dimension greater than three for general g. 

We note how this construction of observables generalizes: The state can be taken to 
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be (j) + = i YJi=\ e ( - < 8 > e, and the observables can be constructed as A x = \ x ■ X and 
By = \ y - X, where X is a vector of matrices X, generalizing Pauli matrices in some 
sense (they anticommute, i.e. XjXj + XjX, = 0 for i 7 ^ j). 


2.5 Bell inequalities allow to lower bound the Hilbert space 
dimension 


In the previous example we chose a to be a square matrix, namely a = 1 121 . We 
will now illustrate the role of the dimension of the measurement directions d' with 
an example of a trivial Bell inequality, where d! = 1 suffices to obtain the Tsirelson 
bound. For this example the coefficients are g = 1 ( 2 \ An obvious singular value de¬ 
composition of this identity matrix is to choose V = S = W = 1^ 2) . Just as before we 
can say that a = is a solution to Eqs. (17 1 and (18 1 , thus the bound is achievable 
with d' = 2. But we can also choose a = (1, l) r , which also solves the system of 
equations. In this case the measurement directions are one-dimensional id' = 1), in 
fact they are all equal to 1. Then the expectation value given by the scalar product 
of the measurement directions reduces to the “classical” expectation value of deter¬ 
ministic local and realistic theories given in Eq. (|4]). Both quantum theory and local 
realistic theories can achieve the maximal value of two. This inequality is therefore 
unable to show a contradiction between quantum theory and locality, realism and 
freedom of choice. You might have expected this, since the matrix of coefficients 
does not even contain a negative coefficient, which implies that the maximum value 
is achieved if all outcomes are +1. 

Let us discuss a more interesting example. It is a special instance of the family of 
Bell inequalities discussed by Vertesi and Pal in l28l . You can also find the following 
analysis for the whole family in the supplemental material of 11271 . The coefficients 
are 


8 = 


/ 

1 

1 

1 

1 \ 


1 

1 

1 

1 


1 

1 

1 

1 


1 

1 

1 

1 


1 

1 

-1 

1 


1 

1 

-1 

1 


1 

1 

-1 

1 

V 

1 

-1 

-1 

1/ 


( 22 ) 


Please note, that the columns of g are orthogonal, thus it is easy to find a truncated 
singular value decomposition of g: We can choose V^ ^ =g. = 2y/2t( 4 ^ 

and W = 1^ 4) . One can easily check, that a = -s/211^ is a solution for the ( d x d')- 
matrix a of Eqs. ( |T7j i and ( fi~8j ), so the maximal quantum value of 16 (see Eq. 
is achievable with (d' = 4)-dimensional measurement directions. It turns out, that 
the system of equations is not solvable if we choose d' = 3, i.e. a to be a (4 x 3)- 
dimensional matrix. This has some very interesting physical implications. Since 
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[d! = 3)-dimensional measurement directions do not suffice to obtain the maximal 
value of the Bell inequality, we can conclude from a measured value of Q = 16, that 
our measurement directions were at least four-dimensional. Of course one will never 
measure this value perfectly in experiment, so what one has to do in practice is to 
calculate the maximum of the Bell inequality over all three-dimensional measure¬ 
ment directions (this is analog to the calculation of the classical bound described 
above). If we call this value T$, then any value between Tj and 16 witnesses the 
dimension of the measurement directions to be at least four (see Figure[3]). 

For spin-4 particles, there are three orthogonal measurement directions (orienta- 


i-H 

II 

E-7 

II 



Fig. 3 Depending on the dimension d' of the measurement directions different values Tji are max¬ 
imal for the Bell inequality given by coefficients in Eq. ( |22[ ). An experimentally obtained value Q 
of the Bell inequality inside the shaded area witnesses, that the produced quantum system had a 
greater Hilbert space dimension than qubits (see text). The values are taken from (28). 


tions of the Stern-Gerlach-apparatus), i.e. x-, y- and ^-direction, corresponding to 
the three Pauli matrices (see Eq. (|6]i) and not more. This holds for all quantum sys¬ 
tems with two-dimensional Hilbert space (qubits). Thus if in some Bell experiment 
the value of the Vertesi-Pal-inequality given by the coefficients in Eq. 22 is found to 
be 16 (or larger than Tj), one can conclude that the produced and measured systems 
were no qubits. In particular they were not single spin-1 particles. Please note, that 
this argument is independent of the physical implementation of the source and the 
measurement apparatuses. For this reason the concept is often called device inde¬ 
pendent dimension witness. 


2.6 Satisfiability of the bound can be understood geometrically 


With r = vff' Eq. (17 1 can be written as r T aa T r = 1. This quadratic form defines 
an ellipsoid with semi-axes —j=, —j=—j= where Ai, A 2 ,..., A^ are the eigenval- 

v ^1 V 2-2 vAf 

ues of aa T . Analogously the vectors r' = \Jwffl lie on the same ellipsoid (see 
Eq.fTS}). 

We therefore state, that the singular value bound is obtainable if and only if the vec- 
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Fig. 4 The singular value 
bound is achievable if and 
only if the vectors V„ and 


■J jj^WyP lie on the surface 
of an ellipsoid. These vectors 
and the ellipsoid (here a 
circle) are shown for the 
CHSH inequality. 


Fig. 5 The vectors V, 
/m7. 


(d) 


and 


V Mi ,, y* ) of 8 = 1<2) lie on 

the dotted ellipse. Increasing 
the larger semi-axis while 
keeping the vectors on the 
ellipse leads to the solid 
(degenerate) ellipse in the 
limit. Infinite semi-axes of 
the ellipsoid imply, that lower 
dimensional measurement 
directions (here d' = 1) suf¬ 
fice to achieve the Tsirelson 
bound. 


wi 




tors vjf ! and lie on an ellipsoid. As we mentioned before, in many cases 

(e.g. from the literature), a can be chosen to be proportional to the identity matrix. 
Thus in these cases the vectors lie on a d -dimensional sphere, i.e. for d = 2 they are 
on a circle, which is shown for the CHSH inequality ifTTj in Fig. [4] 

If a is not square or not full rank (i.e. at least one eigenvalue of a is zero), then 
at least one of the eigenvalues of aa T is zero, too. We define the corresponding 
semi-axis to be infinite. 

The measurement directions lie in the image of the linear transformation associated 
with a. Thus the dimension of the measurement directions cannot be larger than the 
rank of a. For g = 1 we show the degenerate ellipsoid with one infinite semi-axis 
corresponding to the solution a = (1, l) r (see above) in Fig. [ 5 ] 


2.7 Changing g without changing the Tsirelson bound 

The parts of the SVD of g which do not correspond to the maximal singular value 
of (i.e. the non-shaded areas in Figure [2} did not appear in our discussion of the 
Tsirelson bound. Therefore any changes of these singular vectors in V and W and 
singular values in S will not affect our analysis. The last is, of course, only true 
as long as these new singular values do not become bigger than the (previously) 
maximal singular value. While this changes the matrix g, i.e. leads to a new Bell 
inequality, the quantum bound remains obtainable and its value remains the same. 
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Fig. 6 Alice and Bob share pairs of particles in a spin-entangled state \f/ and want to violate a Bell 
inequality. They each can measure the spin of their particle along transversal axes with different 
angle relative to the table’s up. Unfortunately they were not able to agree on what “up” means, 
yet, and their local coordinate systems are twisted by a relative angle <p. The text explains that one 
possibility is to measure using (local) angles yi = 45°, Ji = —45°, 73 = 0°, 74 = 90° at Alice’s site 
and 5i = 0° and 82 = 90° at Bob’s site and “rotate” the Bell inequality. 


From the geometric picture we immediately understand, that rotations of the vectors 
vjf ! and which keep them on the ellipsoid (see Figs, j^jandjsj) also do not 

change the value and satisfiability of the singular value bound. 

We give an example to illustrate that the measurement directions can be rotated 
without affecting the singular value bound and its tightness. Consider the CHSH 
test described above, but now Alice and Bob did not agree on a common coordinate 
system before performing the experiment, see Figure[6] Let us assume for simplicity 
that their local coordinate systems are only rotated relative to each other by an angle 
(p around their common y-axis. This angle <p is unknown to Alice and Bob at the time 
of collecting the measurement data. The quantum state is still t(1,0.0.1 ) 7 , 
independent of (p. 

Let us analyze the effect of the relative rotation on the violation of the CHSH in¬ 
equality. The first idea might be to measure the observables of Section [2~4] in the lo¬ 
cal basis and insert the estimated expectation values into the CHSH inequality. For 
a relative angle <p = 0° these observables are optimal, but for an angle of <p = —45° 
Alice and Bob measure in the same direction and their data will not violate the 
CHSH inequality. From the previous considerations we know that it is also possible 
to “rotate” the Bell inequality such that the actually performed measurements are 
optimal for that inequality. This is can be done by applying a rotation matrix to the 
matrix W. However, twisting the original CHSH inequality by 45° gives \/21 (up 
to relabeling of the measurement settings), see Figures [4] and [5] And as it is shown 
in Figure [5] all but one semiaxis of the ellipse associated with a can be chosen to 
be infinite, which is equivalent to the fact that the classical bound and the quantum 
bound coincide. This implies that the inequality given by coefficients g = \/2 11 can¬ 
not be violated. 
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T 


Violation 


1.10 


1.15 


MAMMA / 


1.05 


Angle cp 


0 


7T 


Fig. 7 The ratio of the maximal qauntum and classical value, the violation, is plotted for the Bell 
inequality given by the coefficients of Eq. (23} as a function of the relative rotation of the two 
laboratories <p. 


The trick is to include more measurement directions. If the measurement directions 
of Alice already uniquely define the ellipsoid associated with a, then the rotation of 
the measurement directions of Bob does not change the fact that the Bell inequality 
can be violated. One obvious possibility to achieve this is to add all settings of Bob 
to Alice. We do this for the CHSH inequality (see Eq. (fl6|) and get 



(23) 


V° 1 


If we call the different measurement angles j\ . , 73,74 at Alice’s site and 8 \. 82 at 

Bob’s site we have for a = y/2t that 71 = 45°, 72 = -45°, 73 = 0°, 74 = 90°, = 0° 

and 82 = 90° are optimal measurement settings. The quantum value T = 4 of this 
inequality does not depend on (p , but the classical bound B does. Figure [7] shows the 
violation of the Bell inequality depending on the relative rotation (p. As expected it 
is always strictly larger than one. The maximal violation of 4 — 2y/2 can be obtained 
for <p =kj, where k is an integer number. We remark that if Alice and Bob even do 
not agree on a common coordinate system for the analysis of the data, they still can 
maximize the violation over the angle (p. 

A similar analysis for a general rotation in three dimensions given by three Euler 
angles was done in l29l . Different approaches to Bell inequalities without a com¬ 
mon coordinate system have been described in the literature. We want to mention 
the following strategy. Each party measures along random but orthogonal measure¬ 
ment directions. Afterwards the violation of the CHSH inequality is calculated for 
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all combinations of pairs of measured settings of Alice and Bob. The result is simi¬ 
lar to the one in this section: if the parties measure along more than two directions, 
then one can find a Bell inequality that is violated with certainty l30l . 

A deeper understanding of the correlations between measurements on separated sys¬ 
tems possible according to quantum theory, including the maximal value of Bell in¬ 
equalities, is an aim of ongoing research in the field of quantum information theory. 
In this essay we saw how more measurement settings and higher-dimensional quan¬ 
tum systems can lead to stronger violations of Bell inequalities, e.g. in the context 
of device-independent dimension witnesses or Bell experiments without a shared 
reference frame. The insights gained from these simple examples may help to find 
Bell inequalities well suited for different situations and applications. 
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Appendix 

3 Tsirelson’s theorem carries the Tsirelson bound to Linear 
Algebra 


We now sketch the derivation of Eq. GD following 11271 . It is strongly based on a 
theorem by Boris Tsirelson EQ . It links the expectation values of quantum measure¬ 
ments to scalar products of real vectors. While the full theorem shows equivalence 
of five different ways of expressing the expectation value, we will repeat two of 
them here. 

Remember that in the formalism of quantum theory observables are hermitean op¬ 
erators, i.e. they equal their complex conjugated transpose. And quantum states can 
be described by density matrices, which are convex mixtures of projectors onto pure 
quantum states, with the weights being the probability to find the system in the cor¬ 
responding pure state. This implies that the density matrix is positive and has trace 
one. 

Consider two fixed sets of observables with eigenvalues in [—1,1], {A i ,A 2 , } 

and and a quantum state given in terms of its density matrix p. 

Then the expectation value of the joint measurement of A a and B y , A x 0 B y , is 
E (x.y) = trf/fv C B y p ) according to quantum theory. Tsirelson’s theorem states, that 
there exist real M\ +M? dimensional unit vectors {vi, V 2 , Vm,} and {wj, W 2 , 
such that all expectation values can be expressed as E(x,y ) = v A • w v . This is the di¬ 
rection we need, because it allows us to replace the expectation value in Eq. 0 
by the scalar product of some real vectors. Tsirelson also proved the converse 
direction: given the vectors Vi,V 2 and w|. vv ’2 ..... w*f 2 there exist observ¬ 
ables A | .Aj. and B\. /A -... ,Bm 2 and a state p such that the expectation value 

E(x,y) = tr (A x ®B y p) equals the scalar product v A ■ w v . 

After application of Tsirelson’s theorem Eq. 0, i.e. Y<x,y8x,yE(x,y), takes the form 


M\ M2 M\ M2 M\ +M2 M\ +M2 

E E 8x,y E V *.,W yJ = E E E E v x,i8x,ySij w yJ 

x=ly=l i=l *=ly=l i—l j=l 




(24) 


Here we expressed the scalar product as a matrix product using the M\ + M 2 di¬ 
mensional identity matrix p( M i+ M 2 ) and defined the vectors v and w, which arise 
if one concatenates all v A and w v , respectively. For the decomposition given in 
Eq. (|16jl with a = l( 2 ), for example, Vi = ( 775 ’ TTo) 7 ^ anc ^ V2 = anc * 


k \fi' \fi‘ 


thus v = — ~^) T ■ F rom Eq. ( |24j i we see, that the maximal quantum 

value of the Bell inequality is given by the maximal singular value (the maximal 
stretching factor) of g(g>l( Ml+A E) times the length of the vectors v and w. The ma¬ 
trix g® has the same singular values as g, except that each of them appears 

M 1 +M 2 times. Because the v A and w v constituting v and w are all unit vectors, the 
length of v is y/M[ and the length of w is y/Mi. Putting these factors together we 
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arrive at Eq. {H}. 


4 MATLAB snippets 

function [ T ] = singularvaluebound( g ) 
%SINGULARVALUEBOUND Calculates the SV-bound of g 
% the returned value is a Ts ire Is on bound for 
% the CHSH—type inequality given by g 
T=sqrt ( numel ( g )) * norm(g ); 
end 

function [ a ] = alphamatrix( g ) 

%ALPHAMATRIX Links SVD to measurement directions 
% see PRL 111, 240404 (2013) 

[Ml M2]= size (g); 

[V S W] = svd (g ); 

acc = lE—4; % adjust to numerical precision 
d=sum( diag (S)>=S(1, 1) — acc ) ; 

% the vectors to be normalized by alpha: 

A= [V(1 :M1, 1 : d ); sqrt (M2/M1)*W( 1 :M2,1 : d )] ; 

Q=( A*A ’). ~ 2 ; 

c = pinv (Q)* ones (M1+M2, 1); 
i f sumf abs (Q*c—ones (M1+M2, l))>acc) 

error! ’ alphamatrix : nosol ’ , ’No u solution^alpha^ found. ’); 

else 

X=A’* diag(c)*A; 
if eigs(X,l , ’ sm ’ )<0 

error( ’ alphamatrix : norealsol ’ , 

’No^real ^solution^alpha^found. ’ ); 

end 

end 

a=X~0.5; 

end 

function [ T ] = tsirelsonbound! g ) 

%TSIRELSONBOUND Calculates the Tsirelson bound for g 
% Uses the semidefinite programm described by 
% Stephanie Wehner in PRA 73, 022110 (2006). 

[Ml M2]= size (g); 

W=[zeros (Ml,Ml) g; g’ zeros (M2,M2)]; 

G=sdp var (M1+M2.M1+M2); 
obj=trace(G*W)/2; 
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F=set (G>0); 
for i = 1 :M1+M2 

F=F+set(G( i,i) == 1); 

end 

solvesdp(F,obj ,sdpsettings(’verbose’ ,0)); 
T=—double ( obj ); 

end 



